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ABSTRACT 

Simple, one-zone models for inhomogeneous chemical evolution of the Galactic 
halo are used to predict the number fraction of zero-metallicity, Population III stars, 
which currently is empirically estimated at < 4 x 10“^. These analytic models min¬ 
imize the number of free parameters, highlighting the most fundamental constraints 
on halo evolution. There are disagreements of at least an order of magnitude between 
observations and predictions in limiting cases for both homogeneous Simple Model and 
Simple Inhomogeneous Model (SIM). Hence, this demonstrates a quantitative, unam¬ 
biguous discrepancy in the observed and expected fraction of Population III stars. 
We explore how the metallicity distribution of the parent enrichment events /(zq) 
drives the SIM and predictions for the Population III fraction. The SIM shows that 
the previously-identified “high halo” and “low halo” populations are consistent with a 
continuous evolutionary progression, and therefore may not necessarily be physically 
distinct populations. Possible evolutionary scenarios for halo evolution are discussed 
within the SIM’s simplistic one-zone paradigm. 

The values of zq depend strongly on metal dispersal processes, thus we investigate 
interstellar mixing and mass transport, for the first time explicitly incorporating this 
into a semi-analytic chemical evolution model. Diffusion is found to be inefficient for all 
phases, including the hot phase, of the interstellar medium (ISM): relevant diffusion 
fengths are 2-4 orders of magnitude smaller than corresponding length scales for 
turbulent mixing. Rough relations for dispersal processes are given for multiphase ISM. 
These suggest that the expected low-metallicity threshold above zero is consistent with 
the currently observed limit. 

Key words: stars: abundances — ISM: abundances — Galaxy: formation — Galaxy: 
halo — galaxies: evolution — early universe 


1 INTRODUCTION 

The search for zero-metallicity, Population III stars provides 
one of the most fundamental empirical constraints on the 
assembly of galaxies and the process of chemical evolution. 
In spite of decades of methodical searches, to date no bona 
fide Population III stars, i.e., having zero metallicity, have 
been identified (e.g.. Beers 1999). This lack is often cited as 
problematic (e.g.. Bond 1981; Cayrel 1996; Larson 1998): a 
first generation of stars must have existed, and to date, the 
stellar initial mass function (IMF) empirically has appeared 
rather robustly universal. We would therefore naively expect 
some surviving Population III stars today. Yet, the on-going 
efforts of Beers et al. (e.g., 1992; 1999) show a current upper 
limit to the fraction of Population III stars in the Galac¬ 
tic halo of < 4 X 10“^ (see §1.1 below). However, there is 
surprisingly little quantitative discussion on the expected 
relative fraction Tm of Population III stars. Bond (1981) 
first emphasized the existence of a discrepancy in the ob¬ 


served fraction of low-metallicity stars for the Galactic halo 
population, in terms of the Simple Model for galactic chem¬ 
ical evolution. A few additional estimates for Fin have been 
made, and these vary by orders of magnitude, ranging from 
Fin ^ 10~^ (Bond 1981) to 10“"^ (Tsujimoto etal. 1999). 
A primary objective of this paper is to further discuss the 
parameters that determine Fm and make better-understood 
predictions. 

Meanwhile, these searches for the most metal-poor stars 
in the Galaxy have more clearly determined their metallicity 
distribution. At present, the lowest-metallicity stars identi¬ 
fied thus far have logarithmic abundance relative to solar 
of [Fe/H] ~ —4. There has been a parallel effort obtaining 
abundance measurements of the lowest-metallicity Lyman 
a absorption-line systems (e.g., Lu etal. 1996; Songaila & 
Cowie 1996). These are currently in the range of [Fe/H] 
~ —3 (Lu etal. 1996; Prochaska & Wolfe 2000; Ellison 
et al. 2000). Pettini (2000) emphasizes that metals have been 
detected in all Ly a forest systems that have been searched 
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to date, including those with the lowest column densities 
(Ellison et al. 2000; Songaila & Cowie 1996). It therefore 
has been suggested that [Fe/H] in the range —3 to —4 repre¬ 
sents a physical minimum threshold in metallicity. Indeed, it 
is widely thought that the abundances of Galactic halo stars 
in this range exhibit the yields of only a few, even individ¬ 
ual, supernovae (e.g., Audouze & Silk 1995; Ryan etal. 1996; 
Shigeyama & Tsujimoto 1998). 

How many Population III stars can we expect to hnd 
today? Besides these, is there indeed a minimum metal¬ 
licity threshold? If so, is its value of order [Fe/H] ^ —4? 
Predictions for these parameters will necessarily have large 
uncertainties, since they depend on poorly constrained re¬ 
lations between cosmic galaxy evolution, star formation, 
and multiphase interstellar medium. Nevertheless, some in¬ 
sight on these fundamental problems is possible from simple 
arguments based on first principles and dominant effects. 
These are re-examined here, with special emphasis on in¬ 
terstellar element dispersal, and lead to new estimates of 
the expected numbers of Population HI stars and the low- 
metallicity threshold. The results also offer insight on the 
star formation processes in the early Universe. 

1.1 Current empirical limit 

We first estimate the current limit on Fm for the Galactic 
halo. Beers (1999) reports that his group’s survey of metal- 
poor stars in the solar neighborhood (e.g.. Beers etal. 1992; 
hereafter BPS) to date finds at least 373 stars with [Fe/H] 
< —2.5. However, none have been found at zero metallic¬ 
ity, or even at [Fe/H] ^ —4. Their survey is not designed 
as a complete survey of the halo metallicity distribution 
function (MDF), thus we compare with the halo MDF mea¬ 
sured by Garney et al. (1996), who hnd 16% of halo stars to 
have [m/H] < —2.5. Their measurements of [m/H], a loga¬ 
rithmic metal abundance relative to solar, are calibrated to 
correspond closely with [Fe/H] (Garney etal. 1987). Thus, 
assuming that all the BPS stars with [Fe/H] < —2.5 are 
halo members, we hnd that the metal-poor tail examined by 
the BPS survey implies a fraction of zero-metallicity stars 
Fill < 4 X 10"% 


2 PREDICTIONS IN THE HOMOGENEOUS 
LIMIT 

What value of Fm should we expect? The very simplest 
generalization is that, for a present-day mean metallicity zi, 
the number of generations of enrichment events is: 

n = zi/zo , (1) 

where zq is the mean metallicity of these individual enrich¬ 
ment events. The present-day fraction of Population HI stars 
is therefore Fm ~ 1/n. Thus, if we wish to interpret the 
metal-poor empirical thresholds as estimates of zo, we can 
set a lower limit on Fm. Taking zq ~ 10~® Zq therefore 
implies n ~ 10® to attain solar metallicity. Likewise, the 
Galactic halo, having a mean metallicity around 0.1 Zq, re¬ 
quires n ~ 10^, and hence Fm ~ lO”"^. 

However, this overly simplistic analysis does not ac¬ 
count for the consumption of gas to form stars, assuming a 
closed system. We therefore turn next to the Simple Model 


for homogeneous chemical evolution (e.g., Schmidt 1963; 
Pagel & Patchett 1975; Tinsley 1980) to describe the re¬ 
lationship between n and metallicity z. The effect of gas 
consumption naturally implies larger numbers of stars for 
earlier generations, thereby increasing the predicted Fm. 

As Bond (1981) demonstrated, the Simple Model pre¬ 
dicts a hnal metallicity distribution: 

/'(2) = y~^ exp{-z/y) . (2) 

For yield j/ ~ 2 i ~ 0.02, this implies that about 15% of 
Galactic halo stars should have [Fe/H] < —2.5. At the time 
of Bond’s work, none of 129 halo objects, namely stars from 
the Garney et al. (1979) proper-motion survey and halo glob¬ 
ular clusters, had been found to have [Fe/H] < —2.6. Ten 
years later, Ryan & Norris (1991) reported that, within the 
errors, their newer dataset did show agreement between the 
prediction and observations; and another five years later, 
Garney etal. (1996) report 16% of stars in their expanded 
survey have [Fe/H] < —2.5, a value in remarkable agreement 
with the prediction of equation 

Is Bond’s discrepancy resolved? Let us now consider the 
number of stars expected below the current observational 
limit of [Fe/H] < —4. Equation predicts a fraction 5 x 10~® 
of halo stars to have metallicities below this value. This is 
about an order of magnitude greater than the observed limit 
found above: about 12 stars should have been found by the 
BPS survey, with assumptions in §1.1. Thus it appears that 
Bond’s discrepancy has simply shifted to lower metallicity. 
While our search techniques are greatly improved, the saga 
of Bond’s original discrepancy is an important caveat against 
prematurely concluding that the discrepancy exists. 

Therefore, the limit on observed zero-metallicity, i.e.. 
Population HI, stars themselves offers perhaps a more pow¬ 
erful constraint on the early chemical evolution and star 
formation conditions. However, it is difficult to use the Sim¬ 
ple Model to robustly predict the actual fraction of zero- 
metallicity stars: equation shows that for low metallicities, 
the predicted number of stars per unit 2 : is roughly con¬ 
stant. Because complete homogeneity is assumed, the Sim¬ 
ple Model cannot accommodate the stochastic effects that 
dominate the first generations of contamination. We there¬ 
fore require an inhomogeneous model for chemical evolution 
that can account for the quantization and stochastic nature 
of the earliest enrichment events. 


3 INHOMOGENEOUS LIMIT: NO MIXING 

As mentioned above, the element abundances of the most 
metal-poor stars are generally interpreted as reflecting the 
yields of the hrst generations of stars. The impetus is par¬ 
ticularly provided by (McWilliam etal. 1995): the marked 
change in abundance patterns for [Fe/H]< —2.4, as is ex¬ 
pected when these properties are determined stochastically 
rather than statistically; and the large dispersions in these 
properties, that are qualitatively consistent with stochas¬ 
tic supernova (SN) yields. As emphasized by McWilliam 
etal. (1995) and Ryan etal. (1996), these abundance dis¬ 
persions now force us to abandon the Simple Model and 
consider inhomogeneous chemical evolution for the lowest- 
metallicity populations. 

In the limit of no ISM mixing, inhomogeneous chemi- 
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cal evolution can be described in the simplest case as over¬ 
lapping regions of contamination. Oey (2000; hereafter Pa¬ 
per I) presents such a model. This model is essentially an 
inhomogeneous variation of the Simple Model, and offers an 
opposing limit for the case of no homogenization. We will 
refer to this as the Simple Inhomogeneous Model (SIM), 
and briefly review some of the basics here. We note that 
Argast et al. (2000) present a similar model that is numer¬ 
ically calculated; their study differs in that they consider 
only individual SNRs, while Paper I considers contamina¬ 
tion from multi-SN superbubbles. In addition, the analytic 
development in Paper I more clearly shows the relevant pa- 
rameterizations that also apply to the numerical results. 

For the metallicity distribution function of long-lived 
stars: 

n n 

n ^^' 

j=l k—j 

where Nj{z) is the MDF for stars born in areas resulting 
from j overlapping polluted regions, is the probability 
of any point having j such overlapping regions after k gener¬ 
ations, and Dk = 1 — kS is the fraction of gas remaining after 
k generations. This assumes that S, the fraction of original 
gas consumed by each generation, is constant; thus it is re¬ 
lated to the present-day gas fraction /ri by = 1 — n5. For 
any given generation k, Pj is given by the binomial distri¬ 
bution: 

= (") Q" (i-O)"-" , i<i<n , (4) 

where Q is the contamination filling factor for each genera¬ 
tion, which is assumed constant. The corresponding proba¬ 
bility of Hnding zero-metallicity gas is: 

Po = (1 - QT ■ (5) 

This model assumes instantaneous recycling, namely, that 
the products of one generation immediately contaminate the 
following. Figure ^ shows a schematic representation of this 
model. 

As discussed more fully in Paper I, the present-day 
mean metallicity for the SIM is characterized by nQ, the 
mean of the binomial distribution. We assume that the 
metallicities for each generation of component enrichment 
events are drawn from a fixed parent distribution f{zo); 
these enrichment events are the units that build up the 
metal abundance. Here, they represent contamination by 
core-collapse SNe originating in OB associations and their 
superbubbles (see §3.2), although the model can accommo¬ 
date f{zo) generated by alternate mechanisms. 

For the Galactic halo. Paper I found that the SIM is in 
good agreement with the observed MDF of long-lived stars 
for a model with nQ ~ 6.4, range logZmin to logZmax of 
-3.7 to -2.0, and present-day gas fraction /ii = 0. Here, 
we further examine kinematic subsets stars isolated by Car¬ 
ney etal. (1996; hereafter CLLA96): their Figures 5, 6, 
and 7 show MDFs for population subsets dehned accord¬ 
ing to kinematic criteria based on, respectively, orbital ve¬ 
locity V, the quadratic radial and perpendicular velocity 

-|- and eccentricity e. These yield very similar 

sequences of MDFs, so we arbitrarily select the sequence in 
eccentricity (Figure 7 of CLLA96) for comparison with the 


SIM models. Figure shows the six subsets of data (dot- 
dashed histograms), with the different ranges of e, increas¬ 
ing in steps of roughly 0.15. The MDFs are converted from 
[Fe/H] to [0/H] according to the relation given by Pagel 


(1989): 




[0/H] 1 

= 0.5 [Fe/H], 

= [Fe/H] -b 0.6, 

[Fe/H] > -1.2 

[Fe/H] < -1.2 

(6) 


Overplotted in solid lines are SIM models with the same 
parameters given above, but with nQ increasing in steps 
of 8, as indicated. These particular models take Q = 0.8, 
varying n in multiples of 10. 

Figure ^ demonstrates that this simple evolutionary se¬ 
quence shows important correspondences with the simple 
sequence in eccentricity isolated by CLLA96: (1) The quali¬ 
tative shape of the SIM MDFs, especially for the hrst three 
subsets in Figures |^ — c, are in reasonable agreement with 
the data. (2) The relative rate in the shape change agrees 
well between the models and data. (3) The models and data 
agree reasonably well quantitatively; in particular, the slow¬ 
ing rate of evolution in the metallicities agree well between 
the models and data. Note that there are important limita¬ 
tions for this SIM because of the particular range for f{zo), 
which is discussed in §§5 - 6. However, as shown in those 
sections, f{zo) and the corresponding SIM can be adjusted. 

The point here is that these fundamental correspon¬ 
dences strongly suggest that the stellar subsets represent 
an evolutionary sequence, and that inhomogeneous evolu¬ 
tion is necessary to interpret the lowest-metallicity sub¬ 
sets. CLLA96 suggest that these subsets represent a transi¬ 
tion from “high halo” to “low halo” populations, such that 
the lower-eccentricity sets represent a disk or “proto-disk” 
population. However, the quantitatively simple progression 
of these eccentricity-dehned subsets as an evolutionary se¬ 
quence could also suggest that these stars are all members of 
a monolithic halo that formed slowly relative to the contam¬ 
ination timescale. Note that the simplistic bimodal formula¬ 
tion of equationjq may exaggerate the lack of low-metallicity 
stars in Figures®—/; inspection of Figures 5 - 7 of CLLA96 
shows a more pronounced low-metallicity tail in [Fe/H]. But 
if this apparent “G-dwarf Problem” in these subset MDFs is 
real, it could be a manifestation of the disk G-dwarf Prob¬ 
lem, as might be expected if these populations are indeed 
related to the disk. 

Regardless of the origin of these kinematic subsets, 
their identification as an evolutionary sequence is strongly 
supportive of the SIM pattern for chemical evolution. The 
shape of the MDF at later evolutionary stages seen in Fig¬ 
ures 1^ — / is essentially characteristic of homogeneous Sim¬ 
ple Model. However, at the earliest stages, the SIM MDF 
shows a distinctly different shape that is influenced by the 
parent f{zo) (see Paper I and §§5 - 6 below). This shape 
must then transition to that of the Simple Model. That this 
is apparently seen in the kinematic subsets therefore strongly 
suggests the SIM behavior and again supports the different 
SIM interpretation of the halo MDF, in contrast to that of 
the Simple Model (see Paper I and §6.3 below). Thus we see 
that the SIM is especially relevant to the lowest-metallicity 
regimes and hence it is appropriate for investigating the fol¬ 
lowing problems. Note that we therefore consider only core- 
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Figure 1. Schematic illustration of the Simple Inhomogeneous Model described by Oey (2000), for n = 2. Each generation k randomly 
occupies filling factor Q. Here, A: = f is shown by white regions and fc = 2 by gray regions. Overlapping regions for j = 2 are shown in 
black, while j = I corresponds to both white and gray regions. 


collapse SN products with the instantaneous recycling ap¬ 
proximation. 


3.1 The fraction of Population III stars 

This inhomogeneous model provides straightforward predic¬ 
tions for the fraction of Population III stars, since the chem¬ 
ical evolution is characterized by the contamination filling 
factor Q and number of overlapping generations n. Note 
that in the case of the Simple Model, only the very first 
generation of stars corresponds to Population III; whereas 
for the SIM, all regions with j = 1 are included for Pop¬ 
ulation III stars. Hence, additional generations contribute, 
provided that some of the star formation takes place in any 
remaining pristine gas. Thus, ifiii is given by: 

n inn 

Tin = Dk-iPi,k /EE Dk-iPj,k , (7) 

k=l J —1 k=j 

where the numerator describes the probability of star for¬ 
mation occurring in primordial gas over each generation k, 
and the denominator corresponds to star formation at all 
degrees of contamination j over all generations (cf. equa¬ 
tion ^). In the n = 2 example in Figure Population III 
corresponds to star formation in the white and gray regions. 
The predicted value of Tm for a given MDF depends only 
on the product nQ and is highly robust 5% variation) to 
different component n and Q. 

For the halo model in Paper I with nQ = 6.4, numerical 
evaluation of equation ^ gives a corresponding prediction of 
Tin = 0.28 for these parameters. Paper I also shows excel¬ 
lent agreement for a SIM of the Galactic bulge, using the 
same parameters as for the halo, but with a more evolved 
nQ ~ 44; this yields Tm = 0.05. These values for Tm are 


again startlingly large and emphasize a blatant discrepancy 
with observations. 

3.2 The low-metallicity threshold 

In the SIM, the present-day mean metallicity zi, deter¬ 
mined by nQ, depends on f{zo), which is characterized by 
its mean, zq. Paper I estimated /(zo) oc 2 “^, within the 
limits^ log 2min and log Zmax of -3.7 and -2.0, respectively, 
yielding a mean log 20 = —3.1. These values were based 
on the simplistic assumption that metal yields of SNe and 
their progenitors are uniformly dispersed within the associ¬ 
ated superbubble volumes, with a power-law distribution of 
in the number of SN progenitors N* per association 
(Oey & Clarke 1998). This is qualitatively consistent with 
observations by Hughes etal. (1998) showing that smaller 
supernova remnants (SNRs) show higher metallicities than 
larger ones. The minimum 2min corresponds to metal dilu¬ 
tion into the largest superbubbles, taken to have radii 1300 
pc and gas density 0.5 cm“®. However, 2inin = —3.7 cor¬ 
responds to [Fe/H] ~ —2.6, which is much larger than the 
observed minimum threshold. This is therefore a significant 
discrepancy between the observations and the SIM having 
these parameters. 

4 INTERSTELLAR DISPERSAL 

The above SIM model represents a limit in the case of no 
large-scale mixing beyond the superbubble radii. The as- 

^ Note in Paper I the values of logz^jj^ and log2niax were inad¬ 
vertently given for niy = 1 Mq instead of the stated rriy = 10 Mq; 
for the quoted value of my = 10 Mq, log2n,in and log2max are 
—3.7 and —2.0, respectively. 


[0/H] N,,, [0/H] [0/H] 
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Figure 2. Kinematic subsets of stars with different ranges in eccentricity e, from Carney etal. (1996; dot-dashed lines), overplotted with 
Simple Inhomogeneous Models having nQ as shown (solid lines). See text for additional SIM parameters. 
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sumption is that metals are uniformly distributed within the 
volumes of the hot superbubbles and simply cool in situ, re¬ 
maining in place over timescale Tn. The values of zo would be 
reduced by dilution in the case of mixing, thereby possibly 
matching the observed lowest metallicities. Since the mean 
Zo is also likely to decrease, a corresponding increase in nQ 
would be required to attain a given present-day metallic- 
ity. This would therefore also reduce the present fraction of 
Population III stars. What are realistic lower limits for the 
values of Zmin and zo? The metallicities f{zo) of typical en¬ 
richment units depend on how far the stellar products of one 
generation are spatially diluted within the inter-generation 
timescale or duty cycle r„. 

We now consider the interstellar mixing process in ordi¬ 
nary, multiphase ISM. Although this is difficult to constrain, 
some hrst attempts have been made by, e.g., Bateman & Lar¬ 
son (1993); Roy & Kunth (1995); and Tenorio-Tagle (1996). 
Here, we examine the mixing process more closely. We con¬ 
sider that there are essentially two mass transport processes: 
diffusion and turbulent mixing. Roy & Kunth (1995) name 
several additional phenomena such as galactic rotation and 
shear, expanding superbubbles, and hydrodynamic instabil¬ 
ities; however, these would all contribute generally to turbu¬ 
lent mixing, and indeed are considered to be the sources for 
interstellar turbulence. Cloud kinematics and collisions have 
also been discussed by Bateman & Larson (1993) and Roy & 
Kunth (1995); however, these will only dominate if and when 
the metal-enriched gas is entrained into cool clouds, and if 
the cloud velocities exceed those of random ISM turbulent 
velocities. Therefore, we will not consider cloud collisions 
here, and the reader is referred to the mentioned studies for 
further discussion of that process. 


4.1 Diffusion 


Diffusion has been discussed in general terms by Tenorio- 
Tagle (1996), who hnds diffusion to be efficient for coronal 
gas typical of the hot ionized medium (HIM), but inefficient 
for cooler gas such as the warm ionized medium (WIM), 
warm neutral medium (WNM) or cold neutral medium 
(CNM). In order to quantitatively compare the efficiency of 
diffusion with turbulent mixing, we evaluate diffusion with 
the Maxwell-Chapman theory, in a more detailed examina¬ 
tion than the approximation used by Tenorio-Tagle. 

Following a point deposition, the concentration, or mass 
fraction of the solute species, evolves under diffusion as, 

where r is distance from the origin point. Mi is the total 
mass of the solute, p 2 is the mass density of the ambient 
held particles, and D is the diffusion coefficient. Following 
the formalism of Woods (1993), the coefficient of mutual 
diffusion for two species of different mass is. 


16 MnoOo ’ ^ ^ 

where T is the thermal temperature; M = mim 2 /(mi -|-m 2 ) 
is the particle reduced mass with mi and m 2 representing 
the masses of the solute and held particles, respectively; 
no = ni + n 2 is the total number density, with ni and n 2 


respresenting the solute and held number densities, respec¬ 
tively; and k is Boltzmann’s constant. The particles are as¬ 
sumed to interact according to an inverse power-law force 
Ki 2 d~'', where d is the distance between particles. 

For ions, the particles repel under the Coulomb force 
with V = 2, for which. 


Ki2 = Zi^2e^/(47r£o) 


( 10 ) 


where Zi and Z 2 are respectively the solute and held par¬ 
ticle charges in units of electronic charge e, and eo is the 
permittivity of free space. We also have, 


Dm = 27 r In A 



( 11 ) 


where A is the ratio of the Debye length to the average 
particle collision impact parameter: 



/ \ZiZ2\e^ 
I A'KeoMg'^ 


( 12 ) 


Te and Ue are the electron temperature and density, respec¬ 
tively; and g is the relative speed between the solute and 
held particles. 

For neutral particles, we can take n = 5, and (Jeans 
1916): 

Ki 2 = kT{v - 1) , (13) 


where a is the distance between particle centers at collision 
for the two species. For atomic H and O interactions, a ~ 
3.15 X 10“® cm. We also have (Woods 1993), 



where Ai(5) = 0.422 and r(|) ~ 1.329. 

From equation ^ we now evaluate the diffusion con¬ 
stants for mutual diffusion of O and H in different inter¬ 
stellar phases, which are given in Table |^, along with the 
assumed temperature, density, and O ionization stage. The 
corresponding r.m.s. diffusion length for a timescale t = 10® 
yr, given by, 

Dms = r2 = QDi2t , (15) 

is also shown. For small concentrations of the solute, D 12 
is insensitive to temperature gradients (Landau & Lifshitz 
1987), so we do not consider a temperature differential be¬ 
tween the two species. Equation ^ shows that D 12 is also 
insensitive to the abundance of O, provided that ni ^ 712 ; 
here we take ni = 10“® 712 . We take Te and Ue to be the 
same as the thermal temperature and density for the H ions. 

The values of D 12 shown in Table are even lower than 
those estimated by Tenorio-Tagle (1996), especially when 
considering the generously low densities assumed here. Most 
significantly, Table shows that even diffusion in the HIM 
is much slower than previously suggested: the scale length 
after 10® yr is only 30 pc, more than two orders of magnitude 
smaller than estimated by Tenorio-Tagle. Thus it appears 
that diffusion is rarely efficient for mass transport. 
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Table 1. Parameters for O diffusion in diffuse ISM 


Phase 

T 

(K) 

n(H) 
(cm ®) 

O ion 

log Di 2 
(cm^ 

' rms 

(pc) 

CNM 

1 X 10^ 

1.0 

O® 

18.25 

0.06 

WNM 

8 X 10® 

0.3 

o° 

20.72 

1.0 

WIM 

1 X 10“* 

0.1 

0+2 

18.04 

0.05 

HIM 

1 X 10® 

0.003 

0+5 

23.64 

30 


^Diffusion length for 1 X 10® yr. 


4.2 Turbulent mixing 

Therefore, it appears that turbulence must strongly domi¬ 
nate the interstellar mixing process in all phases. A mixing 
length scale is given by Bateman & Larson (1993): 

ftrb = ( 2/3 Utrb ^trb , (16) 

where utrb and Ztrb are the characteristic turbulent velocity 
and associated correlation length, respectively. Determin¬ 
ing an appropriate value for Lrb is problematic, since the 
Hi spatial power spectrum in the ISM generally shows clean 
power-law slopes. These have been studied in the Small Mag¬ 
ellanic Cloud (Stanimirovic & Lazarian 2001), Large Mag¬ 
ellanic Cloud (Elmegreen etal. 2001), and the Milky Way 
(Dickey et al. 2001) and found to be broadly consistent with 
Kolmogorov turbulence over the scales considered here. The 
length scale most effective at dispersing a coherent structure 
like a superbubble would be that associated with the object. 
Thus, if we consider Lrb = 50 pc, over the same timescale 
t = 1 X 10® yr, we obtain: 

rtrb = 58.3 (wtrb/kms“^)^^^ pc . (17) 

As a rough first estimate for utrb, we take the soundspeed in 
the given ISM phase, since this represents a characteristic 
velocity. Thus equation shows that for utrb of order 1 — 
100 kms“^, rtrb is 2 - 4 orders of magnitude greater than 
the diffusion lengths estimated in Table 


4.3 The Dispersal of metals 

We now consider the volume over which the metals pro¬ 
duced within superbubbles may be dispersed within the star- 
formation duty cycle r„. Since the various ISM phases have 
different densities and different associated soundspeeds, the 
effectiveness of turbulent mixing depends strongly on the 
phase through which the material is propagating. This in 
turn depends on the cooling process. 

Mac Low & McCray (1988) considered the radiative 
cooling of SN-driven superbubbles, deriving an expression 
for the cooling time trad in terms of the input mechanical 
power Las, ambient number density no and metallicity 

Lad 16 Myr , (18) 

where Las and are in units of 10®® ergs“^ and Zq, re¬ 
spectively. They found that the radiative cooling time is in 
general longer than the superbubble evolutionary timescale, 
for solar metallicities. Since we are considering extremely 
low metallicity conditions, equation demonstrates that 
the first generations of superbubbles should all contain hot, 
coronal gas after the final SNe explode: for ( = 0.1, trad 
already increases by a factor of ~40 from its duration at 

C = 1- 


Heat transport processes are largely analogous to those 
for mass transport, hence the above analysis implies that, 
for hot (10® K) gas, the dominant cooling mechanism is tur¬ 
bulent mixing into the cool ambient medium, rather than 
radiation. As described by Landau & Lifshitz (1987), ther¬ 
mal homogenization by turbulence takes place through the 
mechanical mixing of fluid parcels. Rapid energy dissipation 
by thermal conduction ensures that temperature fluctua¬ 
tions on microscopic scales are always smaller than those 
on larger scales. Thus, the turbulent mixing timescale for 
mass transport should also apply to that for heat transfer. 
However, if the HIM phase dominates the ISM volume, then 
cooling by turbulent mixing may not be efficient (see below). 

Clearly these turbulent processes are extremely diffi¬ 
cult to constrain. We therefore make an attempt here to 
only estimate the order of magnitude scale of these effects. 
As a start, we consider the cooling of material from HIM 
to WIM/WNM (hereafter WM). We can roughly estimate 
the cooling time as the turbulent mixing time for a quantity 
of hot gas mixing with 10 times that quantity of cool gas, 
in number of particles. This ratio implies cooling of the hot 
gas to an equilibrium temperature that is of order 10® K, 
from which the plasma can efficiently cool further by colli- 
sional ionization of He and H (Sutherland & Dopita 1993), 
and metals if any are present. We can estimate the turbulent 
mixing time from equation for hot gas within a superbub¬ 
ble having a final radius R, mixing with cooler ambient gas. 
We take the density ratio between the HIM and WM to be 
ui/u -2 = 0.01, so for for mixing into ambient WM, the rele¬ 
vant turbulent mixing time tcooi is that required for mixing 
the hot superbubble gas with one tenth its same volume of 
ambient WM, thus a total volume with an equivalent radius 
rcooi = 1.03i?. For hot gas mixing into CNM, the equivalent 
i"cooi = 1.0037Z. More generally, we can write, rcooi = eR- 
As argued above, the value of Lrb most effective for mixing 
the superbubbles should be similar to R. This assumption 
is also broadly consistent with the paradigm that superbub¬ 
ble kinematics are a significant driver of ISM turbulence. 
Setting Itrh — R we thus obtain. 


tcool - 2 


Pr 

'^^trb 


(19) 


Again associating utrb with the soundspeed as above, this 
gives a turbulent mixing/cooling time for hot gas: 


tcooi.w 0.10 /trb(pc) Myr (WM) (20) 

tcooi.c 0.94 Lrb(pc) Myr (CNM) , (21) 

for Utrb = 16 kms“® and 1.6 kms“® for cooling into WM 
and CNM, respectively. 

We now consider the dispersal of heavy elements, which 
initially are contained by the hot gas within the parent su¬ 
perbubble. If the breakup of the shell allows immediate ex- 
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posure of the metal-enriched hot gas to the HIM, then tur¬ 
bulent mixing can occur in this highly efficient phase until 
mixing with adjacent WIM and neutral gas cools the hot ma¬ 
terial. Since turbulence dominates both the mass and heat 
transport, we can use limits set by turbulent mixing param¬ 
eters to make a rough first estimate of the dispersal length 
scale for mixing within the HIM. We first consider the frac¬ 
tion of metal-enriched gas that has cooled after an interval 
t: 

where V 2 and (.2 are parameters for WM or CNM, depending 
on the assumed ambient medium. The remainder is still hot 
and can reach a distance rtrb according to equation p^ , with 
4rb = R as before. The total volume ultimately occupied by 
dispersing hot gas is therefore, 

P'^cool 

Vi= (1 - /cool) 47rr?rb drtrh ■ (23) 

Jo 

Integrating this, we find the equivalent length scale 
(3Vi/47r)i/3: 

rmix = i? , (24) 

where subscripts 1 and 2 denote parameters for the dis¬ 
persing and ambient media, respectively. For a density ra¬ 
tio ni/n 2 = 0.01 and vi/v 2 = 10, as we assume for HIM 
{v 2 = 160 kms“^) mixing into WM, equation reduces to 
rmix = 2.6 R for the characteristic mixing distance. 

The extent of mixing within the WM itself is difficult 
to estimate analytically. For a continuous WM, once the hot 
gas has cooled to this phase, individual particles could in 
principle be transported within this phase during long time 
scales. Conversely, cooling in this phase takes place rapidly 
if its global, or local, heating sources are not maintained. 
However, it again seems likely that turbulence determines 
the circumstances that permit individual gas parcels to cool 
from this phase. Considering now the analogous dispersal of 
WM into CNM, we again have ni/n 2 — 0.01 and £2 = 1.03, 
therefore also resulting in equation ^ We adopt the original 
value of R; while the geometric distribution for the cooled 
superbubble gas is greater than R, the cooled material is 
now at the WM density and its equivalent volume is reduced 
by a factor of 0.01 for the adopted densities. Since R is the 
characteristic length, it remains a reasonable length scale for 
equation ^ within the large uncertainties of this approxi¬ 
mation, and also helps offset the likely additional time or 
distance distance needed to encounter CNM. 

Thus, in a first crude estimate, we take the efficiency 
of element dispersal in the WM to be the same as in the 
HIM. Although turbulent mixing itself is less efficient, this 
is counteracted by the likely longer residence time within 
this phase if a continuous WM exists, as it does in local 
star-forming galaxies. 

Meanwhile, mixing also takes place in the CNM. As¬ 
suming mixing of hot superbubble gas with ambient CNM 
starts at tf, the time at which the shell attains its final, 
pressure-confined radius, equation ^ gives a length scale 
for mass transport within the CNM: 


-V2R {Tti tcool,C t/) 


1 1/2 


(25) 


where V 2 is now the turbulent velocity for the CNM. Oey & 
Clarke (1997) give for tf. 

tf = te R/Re , (26) 


where tf. ~ 40 Myr is the duration of SN mechanical power 
for the superbubble, and 


^ 71/2 P2 


-1/2 pl/2 


(27) 


with P 2 as the ambient total pressure. 

It is again essential to bear in mind that these estimates 
are crude approximations, with transport within the WM 
especially uncertain. Full hydrodynamical simulations are 
necessary to more reliably constrain the mixing processes. 
However, it is encouraging that first attempts at numeri¬ 
cal results are broadly consistent with the above estimates 
(de Avillez & Mac Low 2002). In general, the parameters 
adopted above should tend toward upper limits on the dis¬ 
persal length; for example, adopting the soundspeed for utrb 
is generous, and the influence of magnetic fields will gener¬ 
ally inhibit mixing. 

Figure ^ shows the dispersal length rmix as a function 
of superbubble final radius R as given by equations ^ and 
p^ . The HIM/WM dispersal shows a linear relation between 
rmix and R with slope 2.6 (upper set of lines), as determined 
above, and the slope for the CNM limit is ~0.7 (lower set 
of lines). These represent crude limits to the dispersal of 
elements by turbulence. Estimates are shown for star for¬ 
mation duty cycle times of = 500 Myr (dashed lines) and 
200 Myr (solid lines). These values of limit the largest R 
to roughly 500 pc and 200 pc, respectively, as indicated by 
equation Eh; larger objects may exist, but will not be able 
to cool within time Tn. Thus, for characteristic over the 
relevant period of star formation, the superbubbles have a 
corresponding range in R. 

It is apparent that the dispersal length is sensitive to the 
ISM phase balance, with the slope of the relation between 
rmix and R ranging from about 0.7 to 2.6. For normal or low 
star formation rates, this suggests that the dispersal scale is 
a factor of up to 2.6 times the parent superbubble radius, or 
a factor of up to ~20 times the parent superbubble volume: 
i.e., roughly an order of magnitude for this crude analysis. 

Beware that we assume that cooling of the hot super¬ 
bubble gas by turbulent mixing determines the HIM resi¬ 
dence time. Estimates for most nearby galaxies having nor¬ 
mal star-formation rates show low HIM filling factors (e.g., 
Oey etal. 2001). However, if the HIM filling factor is large, 
then the newly-produced metals could be transported in¬ 
definitely. Mixing would take place extremely efficiently in 
a HIM-dominated ISM, justifying instantaneous mixing ap¬ 
proximations in chemical evolution models. On the other 
hand, a HIM-dominated ISM probably implies galactic su¬ 
perwinds, complicating the fate of the stellar products. In¬ 
deed, a blow-out or galactic superwind would return us to 
the conventional problem of losing an unknown amount of 
metals from the system. 


5 INHOMOGENEOUS LIMIT: DISPERSAL BY 
MIXING 

In short, the preceding section gives crude arguments that 
suggest that the dilution of metals beyond the parent su- 
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Figure 3. Estimates for the dispersal length rj^ix from superbubbles having final radii R. The limiting cases correspond to dispersal of 
hot superbubble gas into HIM and/or WM, and CNM as indicated; dashed lines show models with = 500 Myr and solid lines show 
T = 200 Myr. 


perbubble radii is roughly another order of magnitude by 
volume. Note that if the hot, metal-bearing gas is vented 
out of a gaseous galactic disk, that these arguments present 
an upper limit to extended dilution within the disk where 
subsequent stars presumably form. Moreover, merging be¬ 
tween adjacent contaminated regions is greatly enhanced by 
the mixing process, and limits the dilution of metals within 
a galaxy’s star forming volume. 

We now return to the Simple Inhomogeneous Models, 
bearing in mind that these represent a limit for the case of 
no homogenization, i.e., mixing is limited only to within the 
rmix argued above. We can reexamine the values of Zmin and 
zo of individual enrichment events, as well as enrichment by 
the first generation of stars. In §3.2, we had log «min = —3.7 
and logzo ~ —3.1, based on the range of zq estimated in 
Paper I for no interstellar mixing beyond the superbubble 
radii. Since the relationship between Cmix and R is virtually 
linear, we can simply scale the expected values of zq from the 
no-mixing values, thus decreasing them by up to an order 
of magnitude, as estimated above. Note that if the super¬ 
bubble filling factor Q is high enough, mixing from adjacent 
regions will overlap, so this is again a conservative estimate 
in maximizing the decrease in zq. 


5.1 Mixed-dispersal, high-evolution SIM 

Scaling the original range in zq by a factor of 10 gives 
—4.7 < logzo < —3.0 and log2o ~ —4.1. For most chem¬ 
ical evolution models, the reduction in zq implies that the 
number of enrichment generations n must be increased by 
roughly the same factor to achieve a given present mean 
metallicity. With the new range in zq, the SIM in §3 fits 
the Galactic halo MDF with nQ — 64, now predicting 
Fin ~ 3 X 10“^. This is an order of magnitude smaller 
than Fill ~ 3 x 10“^ predicted by the no-mixing SIM in 
§3.1. However, an important consequence of reducing zq here 
is that the shape of the predicted MDF no longer agrees 
as well with the observations as the original model with 


log.zo = —3.1. The dot-dashed histogram in Figure ^ shows 
the observed halo MDF from Carney etal. (1996), using 
their field star sample of 135 stars having retrograde veloc¬ 
ities (V < —220 kms“^). The data are converted to [0/H] 
from [Fe/H] according to equation^ The evolved SIM with 
nQ = 64 is overplotted with the solid curve. We have mod¬ 
ified the conversion from « to [0/H] used in Paper I: where 
previously we simply used [0/H] = log a-I-1.7, we now adopt 
a more exact relation, [O/H] = log[^/(24.32 —128 Z)]-|-3.07. 
The metallicity unit in the SIM is the mass ratio of metals 
to H and He, thus z = Zjil — Z), where Z is the conven¬ 
tional mass fraction of metals. The SIMs assume a default 
SN yield of 10 Mq and present-day gas fraction = 0. 

Note again that this high-evolution case accounts for 
the complete dilution of metals from each OB association, 
but it does not account for mixing between the discrete en¬ 
richment regions, which corresponds to homogenization. As 
discussed in Paper I, such homogenization would cause the 
models simply to tend toward the Simple Model; hence in 
principle, the true distributions should be bounded by the 
Simple Model and the SIMs. Figure^ also shows the Simple 
Model for pi = 0 (dotted line). 


5.2 Mixed-dispersal, low-evolution model 

In §3.2 we saw that the low-metallicity tail for the no-mixing 
model does not extend to metallicities as low as those ob¬ 
served. (Note that although in Figure ^ the CLLA96 data 
are not shown to extend below [O/H] ~ —3, the current ob¬ 
served limit does extend to the equivalent of [O/H] ~ —3.4.) 
We can force the SIM to match the observations by modi¬ 
fying the distribution f{zo), extending only Zmin by a factor 
of 10, so that —4.7 < logzo < —2.0. The upper limit repre¬ 
sents contamination from single SNe, thus individual yields 
and stochastic mixing effects could plausibly maintain a high 
Zmax- Keeping the other input parameters the same, it re¬ 
mains possible to match the observed halo MDF. Owing to 
lack of computational resolution, we cannot run our exact 
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Figure 4. Observed Galactic halo MDF from Carney etal. (1996; dot-dashed line) overplotted with the inhomogeneous models for mixed 
dispersal. The solid histogram shows the low-evolution model (zq = —3.7, nQ = 24), and the solid curve shows the high-evolution, model 
(zq = —4.1, nQ = 64). The Simple Model is also shown, with the dotted line. 


numerical models beyond nQ ~ 30, but the behavior of the 
allowed models shows that nQ now needs to be evolved up 
to nQ ~ 50 to match the data. In lieu of this model using 
the full required range of zo, Figure^ shows a similar model 
with a range —4.4 < log zo < —2.0 and nQ = 24 (solid 
histogram). Comparison of this model with the no-mixing 
model in Figure 2 of Paper I clearly shows that extending 
the lower range of zq will allow an excellent fit to the ob¬ 
served halo MDF, including the low-metallicity threshold. 
This modification may also be applied to the modeled halo 
subsets shown in Figure (again, computational limitations 
prevent our explicit modeling of these subsets). 

The high-metallicity tail of the MDF for this model, as 
for the no-mixing model, reflects the remnant power-law dis¬ 
tribution of /(zq) (Paper I); whereas for the model presented 
in the previous section, the MDF has evolved to sufficiently 
high metallicities that the form of /(zo) is no longer man¬ 
ifest. We therefore refer to the former as “low-evolution” 
models and to the latter as a “high-evolution” model. Al¬ 
though, for the case with extended range in zq, the evolu¬ 
tionary parameter nQ ~ 50 is numerically similar to that in 
the high-evolution model, we still designate this as a low- 
evolution model, since the tail of the parent /(zq) is still 
apparent. For the low-evolution models, equation ^ must be 
evaluated numerically, whereas the high-evolution models 
can be computed with an analytic approximation (Paper I). 

The agreement between the data and the SIM does 
depend on /(zo) ranging up to roughly logzo ~ —2.0, 
which corresponds to [Fe/H] ~ —0.6. This causes the high- 
metallicity tail of /(zo) to remain manifest until the mean 
z evolves to comparable values (see Paper I). Thus the as¬ 
sumed range in zo spans almost three orders of magnitude. 
This range is large, but it is unclear whether it is implausi¬ 
ble (see §6.4 below); the strength of this model is primarily 
that it matches the data. With nQ ~ 50, this model predicts 
Fill ~ 4 X 10 ^. 


6 DISCUSSION 

6.1 The low-metallicity threshold 

The observed low-metallicity threshold, if it is real, provides 
an important constraint on the models by limiting Zmin. As 
we saw above, the model MDFs, evolutionary state, and pre¬ 
dicted Fin are all sensitive to the form and range of /(zo), 
especially Zmin. Fields etal. (2002) analytically demonstrate 
the interrelation between Zmin, the metallicity dispersion, 
and the present-day metallicity. Besides dilution by mixing, 
another model parameter that affects Zmin is the SN yield 
niy. The SIMs assumed niy = 10 M©, which could be over¬ 
estimated by up to an order of magnitude. Reducing niy to 
1 Mq would simply reduce the values of zo, hence zo, by an¬ 
other order of magnitnde. For the mixed-dispersal models, 
the estimated range of the threshold log Zmin thus becomes 
roughly -5.7 to -4.7, corresponding to [Fe/H] ~ —4.6 to - 
3.6. This agrees well with the observed lower limit —4 < 
[Fe/Hj. 

The lowest expected stellar metallicities could conceiv¬ 
ably be below the characteristic Zmin if Zmin varies over time. 
Although these SIMs do not accommodate an evolution of 
/(zo), we briefly consider mixing processes for the earliest 
conditions of star formation. The very first generation of 
stars should take place in WNM since the lack of metals 
in the early Universe maintains temperatures around 500 - 
1000 K at ISM densities (e.g., Barkana & Loeb 2001). As 
subsequent generations of stars begin to generate WIM, the 
production of metals soon enables cooling and the formation 
of a CNM phase. But for the first generation, metal-enriched 
material from the first SNe could mix for an extended period 
of time, since cooling below WNM temperatures is difficult. 
This is also likely to lengthen the timescale for forming the 
next generation of stars, to a value larger than the charac¬ 
teristic Tn at later times. These factors argue for widespread 
dispersal of metals produced by the first generation. On the 
other hand, the character of interstellar turbulence should 
be different from that in the conventional multiphase ISM: 
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since there is no pre-existing stellar mechanical feedback, 
tnrbnlent mixing may be substantially less efficient in dis¬ 
persing the first stellar products. As seen in §4.1, the com¬ 
peting process of diffusion is a much less effective mixing 
agent than turbulence. It is therefore difficult to constrain 
the ISM metallicities resulting from the hrst generation of 
stars. They could be either significantly lower than zo, or, 
in the case of poor mixing, higher than zq. At this stage, 
the value of 2min, taken as a characteristic value over all 
stellar generations, is therefore about the best estimate we 
can make for the low metallicity threshold above primordial. 
As shown above, the rough estimates of Zmin are compatible 
with the currently observed limit. 

We also note that other mechanisms for a low- 
metallicity threshold have been proposed, that do not 
rely on the existence of an intrinsic threshold in the par¬ 
ent enrichment units. For example, Hernandez & Ferrara 
(2001) estimate an effective low-metallicity threshold around 
log Zmin/ Z© ^ —4.3; they obtain this value by estimating 
the number of progenitor sub-halos needed to assemble the 
Galaxy, and extrapolate Zmin as the mean cosmic abundance 
corresponding to individual such sub-halo at the time of 
its formation. Schneider etal. (2002) propose that a sim¬ 
ilar log 2min/ Z© ^ —4 is a physical threshold needed to 
form low-mass (^ 1 M©) stars, with pre-enrichment pro¬ 
vided by a first generation of only supermassive stars. The 
value log Zmin/ Z© of -4.0 corresponds to [Fe/H] of -4.6, thus 
these estimates are quite similar to the value estimated here, 
with large uncertainties admitted by all these studies. Note 
that these alternative approaches avoid addressing the par¬ 
ent f{Z), which is the straightforward origin of Zmin esti¬ 
mated in this work. 

In the above analyses, including our own and those just 
described, Zmin is not necessarily a hard limit, but a charac¬ 
teristic value for the lowest expected metallicities; stochasti¬ 
cally, it is likely that some stars will be found at metallicities 
anywhere between zero and Zmin, depending on their specific 
parent interstellar mixing conditions. However, Zmin should 
correspond to a sharp drop-off in observed MDFs. 

6.2 The fraction of Population III stars 

Besides the Simple Model analysis by Bond (1981), there are 
a couple other predictions in the literature for Fiu, the frac¬ 
tion of zero-metallicity. Population HI stars. Cayrel (1996) 
makes a rough estimate of Fui ~ 10“® for the entire Galaxy. 
This prediction is based on SN contamination within indi¬ 
vidual clouds, and the author emphasizes the uncertainty 
caused by the arbitrarily assumed number of these clouds. 
Tsujimoto etal. (1999) also make a prediction based on a 
detailed model of single SN contamination within clouds, 
with more complex chemical regulation determined by the 
triggering of subsequent star formation in SNR shells. They 
predict Tm ~ 10“® — 10'^ for the Galactic halo. However, 
it is unlikely that the occurrence of star formation is dom¬ 
inated exclusively by triggering in SNR shells, and the ac¬ 
tion of multi-SN superbubbles causes substantial differences 
in ISM dynamics from the assumed individual SNe (e.g., 
Oey & Clarke 1997). Here, we emphasize different parame¬ 
ters, in particular f{Z), thereby presenting a complementary 
perspective on the dominant processes. 

As can be seen above, our Galactic halo predictions for 


Fill are presented in simpler terms, and eliminate most of 
the difficult parameters used by these other authors. Thus, 
the contrast with the empirical limit on Fin should yield 
better-understood constraints on early star formation and 
enrichment. To summarize: the homogeneous. Simple Model 
predicts a fraction 5 x 10“® for stars having [Fe/H] < —4; 
with the no-mixing SIM, we predicted Fin ~ 2 x 10“^; 
considering metal dispersal by mixing, we obtain a high- 
evolution SIM and a low-evolution SIM, both predicting 
Fin ~ 3 — 4 X 10~^. These predictions are all 1 - 3 orders of 
magnitude greater than the observed upper limits. 

The SIM predictions are based on parent enrichment 
events whose metallicity distribution /(zq) oc z“^ is derived 
from dilution of SN products into volumes that are refer¬ 
enced to their SN-driven superbubbles (§3.2 above; Paper I). 
Although, as emphasized above, the mixing process is rather 
uncertain, it is apparent that regions of contamination will 
ultimately begin mixing together and homogenizing, espe¬ 
cially if mixing is even more efficient than the generous es¬ 
timate determined in §4.3. Thus this condition imposes a 
likely limit to the minimum zo driven by dilution of parent 
enrichment events. Once the contaminated regions begin ho¬ 
mogenizing, the system should approach the homogeneous 
limit of the Simple Model. 

We have seen above that Fin depends rather strongly on 
/(zo). Although we used metal dilution within superbubbles 
to derive /(zo) oc z“^, it seems fairly easy to obtain a sim¬ 
ilar power-law distribution in zo via other processes. Thus 
these particular SIMs are fairly robust and independent of 
the aetual origin of the /(zo) distribution. Note also that, 
although we have differentiated the SIMs as “no-mixing” 
and “mixed-dispersal,” the only technical difference between 
them is in the range of zq. As seen in §5, reducing the value 
of Zo implies a reduction in Fin by roughly the same amount. 
Thus if we invoke an overestimate in the default SN yield 
niy by a factor of 10, and take /(zo) for the high-evolution, 
mixed-dispersal model, we obtain minimum estimates for zo 
of —5.7 < logzo < —4.0. This decreases the predicted Fin 
for that model by another factor of 10, to Fin ^ 10“® , 
which is still an order of magnitude larger than the current 
empirical limit. 

Besides /(zo), the SIMs for a given nQ depend only 
on the present-day gas fraction /ri, which was assumed to 
be jzi = 0. Increasing pi will decrease the number of con¬ 
taminating generations required to attain the present-day 
metallicity, hence would only serve to increase Fin. 

Therefore, given that the predictions and observations 
disagree for both limiting cases, the homogeneous Simple 
Model and the SIMs, it seems clear that we can finally 
demonstrate a quantitative, unambiguous discrepancy in the 
observed and expected numbers of Population III stars, based 
on our stochastic, inhomogeneous expansion of the Simple 
Model. 

6.3 Predicted and observed MDFs 

We now consider the predicted and observed MDF’s in their 
entirety. The Simple Model and the no-mixing SIM were dis¬ 
cussed in Paper I. The high-metallicity turnover for the ho¬ 
mogeneous, Simple Model results purely from the consump¬ 
tion of gas to form stars. In contrast, although the MDFs are 
similar, the high-metallicity tail in the low-evolution SIMs 
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reflects the tail of f{zo). On the low-metallicity side, Figure ^ 
shows that the SIMs have distributions similar to the Sim¬ 
ple Model. This is expected, since the SIMs approach the 
Simple Model distribution in this regime (Paper I). How¬ 
ever, unlike the Simple Model, the SIMs have a minimum 
threshold metallicity set by Zmin- As determined in §6.1, the 
predicted low-metallicity threshold is not inconsistent with 
that suggested by the data. 

Since the Simple Model has no such lower limit, it there¬ 
fore predicts too many stars below the observed threshold. 
But aside from this discrepancy, there is good agreement 
between the halo data and Simple Model. This could there¬ 
fore suggest complete homogenization, although paradoxi¬ 
cally, the similar agreement for the low-evolution SIMs like¬ 
wise suggests poor homogenization. As discussed above, the 
mixed-dispersal, low-evolution SIM can produce a minimum 
2 that matches the observed threshold, as well as match the 
observed halo MDF. It therefore best reproduces the data, 
out of the four models. 

Figure ^ shows that the high-evolution SIM does not 
agree with the observed halo MDF, particularly for the 
high-metallicity tail. As demonstrated in Paper I, the high- 
metallicity tail in this model continues to drop off more 
steeply as the system evolves. The disagreement implies that 
this model does not properly account for the dominant ef¬ 
fects. The steep, high-metallicity drop-off characterizes an 
evolved system, as dictated by nQ and zq. Hence, a low- 
evolution model may be a better representation of the sys¬ 
tem. Alternatively, the halo may be better described by fully 
homogenized, rather than inhomogeneous chemical evolu¬ 
tion, thereby implying that the Simple Model best represents 
the system, with the caveat that the lowest metallicities are 
not observed. 


6.4 Halo evolution 

Both the homogeneous and inhomogenous models consid¬ 
ered here are essentially one-zone models, therefore imply¬ 
ing a single coherent formation process for the Galactic halo. 
This is more consistent with the monolithic collapse model 
of Eggen, Lynden-Bell, & Sandage (1962), than the more fa¬ 
vored Searle & Zinn (1978) multi-fragment assembly model. 
Thus, observed discrepancies with predictions found here 
may well be due to a such a fundamental shortcoming in 
the one-zone models. In any case, quantitative characteri¬ 
zation of the comparison and contrast between the simplest 
one-zone models and the data strongly constrain the forma¬ 
tion models. Here we explore some possible implications for 
halo evolution within the one-zone paradigm. 

If we take the observed low-metallicity limit of —4 < 
[Fe/H] at face value, then this constraint implies that 
the Simple Model is not a useful representation of the 
early Galactic halo evolution. We have seen that this low- 
metallicity threshold is, however, compatible with the SIMs 
considered here. The observed lower limit effectively sets 
Zmin for the parent enrichment events, since the observed 
metallicity cannot be lower than Zmin, excepting Popula¬ 
tion HI stars. We saw above that this constraint on f{zo) 
generally yields SIMs for which nQ ~ 50 to 100 to fit the 
observed halo metallicities, i.e., n ~ 50 — 100 generations of 
star formation, for contamination filling factor Q of unity. 


Interestingly, this halo nQ is compatible with solar metal¬ 
licity systems requiring roughly a Hubble time. 

It seems most plausible that early halo chemical evolu¬ 
tion proceeded inhomogeneously, driven by discrete enrich¬ 
ment events in a warm, hard-to-cool ISM with little SN- 
driven turbulence. This is especially suggested by the SIM 
evolutionan' sequence seen in the kinematic subsets shown 
in Figure H. However, we have seen that, above the low- 
metallicity threshold, the Simple Model does match the ob¬ 
served halo MDF as well as the low-evolution SIMs (Fig¬ 
ure ^ . Thus it may be that the early evolution took place 
inhomogeneously, but rapidly became essentially homoge¬ 
neous. As demonstrated in §4.3, mixing takes place rapidly 
if it can take place in a fully developed HIM. For starburst¬ 
like star formation rates, the HIM dominates the ISM, thus 
presumably driving instantaneous mixing, although on the 
other hand, metals are then likely to be lost through galactic 
outflows. 

For halo evolution dominated by inhomogeneous ef¬ 
fects, apparently a large range in zq is required, roughly 
—4.7 < log zo < —2. The lower limit is set by the observed 
threshold, and the upper limit is constrained by the ability 
to reproduce the observed high-metallicity tail. It is unclear 
whether this large range in zo is unrealistic; as shown in §4.3, 
mixing processes that dilute the metals and reduce zo are 
uncertain. Nor do we know whether the mixing processes 
apply to all, or just to some, enrichment events. It does 
not seem obviously implausible that the parent enrichment 
events could span almost three orders of magnitude in zq, 
although the range does seem large. To match the halo data, 
the only real requirement for the low-evolution SIM is that 
f{zo) span this range and follow a decreasing power-law that 
can fit the high-metallicity tail. 

The above interpretations for halo evolution depend on 
the assumption that the observed limit —4 < [Fe/H] rep¬ 
resents a true threshold. Should this not be the case, then 
this would either imply that homogeneous evolution more 
strongly dominates, or that the range in parent zq is even 
larger than implied here. Ultimately, as discussed above, it 
will probably be necessary to invoke a more complex evo¬ 
lutionary process than described by these simple one-zone 
models. In particular, if the halo is assembled from a multi¬ 
tude of individual, smaller star-formation units as suggested 
in the classic Searle & Zinn (1978) scenario, then the co¬ 
herent formation scenario described by these simple models 
presumably would be easily disrupted. 

It is also essential to bear in mind that the discrep¬ 
ancy in the predicted fraction Em of Population HI stars 
poses a profound problem for all these models. A variety of 
solutions has been suggested to explain the absence of ob¬ 
served Population HI stars, for example, a primordial IMF 
composed of high-mass stars, pre-enrichment of the Galactic 
halo system, residence of Population HI stars at large radial 
distances, disguising of Population HI stars by metal accre¬ 
tion or surface convection, etc. The myriad possibilities are 
addressed in the excellent conference proceedings on the first 
stars (Weiss et al. 2000). Many of these suggestions are com¬ 
patible with the simple models discussed above, but again, 
it may be that the discrepancy in Em is symptomatic of a 
more fundamental omission in these simple models of the 
halo evolution. 
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7 CONCLUSION 

We have considered simple, one-zone models for homoge¬ 
neous and inhomogeneous chemical evolution of the Galac¬ 
tic halo to estimate the expected fraction Um of zero- 
metallicity, Population III stars, and to estimate the ex¬ 
pected low-metallicity thresholds above zero metallicity. 
Since the most metal-poor enrichment events depend on 
metal dispersal in the ISM, we also investigated the inter¬ 
stellar mixing and dispersal processes. We then investigated 
the results in terms of the observed halo MDF, offering some 
possible constraints on the halo evolution. Interestingly, the 
observed halo MDF does not itself preclude a monolithic 
halo formation model. Kinematic subsets of stars previously 
identified as “high halo” and “low halo” are broadly consis¬ 
tent with an evolutionary sequence in terms of the Simple 
Inhomogeneous Model. 

Our mass transport analysis finds that diffusion is in¬ 
efficient for metal dispersal since, for the cold and warm 
ISM phases, the relevant diffusion lengths are 2-4 orders of 
magnitude smaller than comparable turbulent mixing length 
scales. This result confirms the earlier approximations by 
Tenorio-Tagle (1996), although with some quantitative dif¬ 
ferences. We also find that diffusion is ineffective even in the 
HIM, for which the relevant diffusion length is still over an 
order of magnitude smaller than the HIM turbulent mixing 
scale. Turbulence also likely dominates cooling processes for 
hot gas. We estimate some rough relations for the dispersal 
of metals in the cold, warm, and hot ISM phases. As with 
diffusion, mixing takes place much more rapidly in hotter 
phases than cooler phases, assuming that turbulent veloci¬ 
ties are linked to the soundspeed. We estimate roughly that 
turbulent mixing dilutes the metallicity of regions originat¬ 
ing as SN-driven superbubbles by up to an order of magni¬ 
tude by volume for ordinary, multiphase ISM. 

This result implies that parent enrichment events may 
have metallicities zo that are up to an order of magnitude 
lower than in the absence of interstellar mixing. Compared 
to the inhomogeneous halo model shown in Paper I, this 
revises the lowest expected values of zo to values that are 
similar to the lowest observed metallicities, —4 < [Fe/H]. 
However, assuming a range in zo of about two orders of 
magnitude, with zq = —4.1, the resulting Simple Inhomoge¬ 
neous Model must be in a highly evolved state to show the 
observed halo metallicities. The corresponding MDF drops 
off too steeply to match the data in the high-metallicity tail. 
This model predicts Fin ~ 3 x 10“^. 

However, as found in Paper I, the observations do match 
a low-evolution inhomogeneous model MDF, especially if the 
range in zq is extended to cover almost three orders of mag¬ 
nitude. This range in values may seem implausibly large, 
but the resulting low-evolution model does best match the 
observed MDF from among the models considered here. The 
required extension in Zmin could be caused by a number of 
effects, for example, an overestimated SN yield and dispersal 
processes. This model predicts Fm ~ 4 x 10~^. 

The homogeneous, Simple Model also matches the ob¬ 
served MDF well, although it naturally predicts no low- 
metallicity threshold, and predicts a fraction 5 x 10“® of 
stars to have [Fe/H] < —4, below the observational limit. 
If the observed threshold is real, it requires a different, pre¬ 
sumably inhomogeneous, model to describe the early halo 


evolution, since large dispersions in metal abundances are 
seen at the lowest metallicities. It may be possible that the 
evolution proceeded more homogeneously at later times. 

Our estimates of Fm with these simple models are all 
at least 1-3 orders of magnitude higher than the empirical 
upper limit of Fm < 4 x 10~^. Given that the observations 
disagree with both limiting cases, the homogeneous Simple 
Model and the Simple Inhomogeneous Models, we therefore 
demonstrate an unambiguous discrepancy in the expected 
number of Population HI stars. This poses a significant prob¬ 
lem for these models, although they remain compatible with 
some of the proposed solutions. Recall that this analysis ex¬ 
amines only the most fundamental processes, and does not 
include parameters like inflow/outffow and time evolution in 
Q, yield, or IMF, any of which could have important effects. 
However, the constraints offered here form the foundation 
for subsequently exploring such effects. 

It is interesting that the current observed lower limit 
of —4 < [Fe/H] is consistent with the lowest-metallicity es¬ 
timates crudely expected for the parent enrichment events. 
Since the lowest metallicities depend on dilution by inter¬ 
stellar dispersal, in principle the metal-poor regime of the 
MDF should also provide constraints on the earliest enrich¬ 
ment processes. Our interpretation of the models depends on 
whether the observed limit is indeed real. Twenty years ago 
we had an apparent discrepancy between the expected and 
observed numbers of stars having [Fe/H] < —2.5; these stars 
have since been found. Today, the Simple Model predicts a 
discrepancy of similar magnitude for [Fe/H] < —4, while the 
inhomogeneous model does predict this metal-poor thresh¬ 
old. Will these additional “missing” stars also be found? It 
is essential that the low-metallicity MDF and threshold, if 
any, be firmly established. 
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